Scalar Boundary Conditions in Lifshitz 

Spacetimes 



Cynthia Keeler 



Department of Physics, University of Michigan, 
Ann Arbor, MI-48109, USA 



keelercOumich . edu 



Abstract 

■ , We investigate the conditions imposable on a scalar field at the boundary of the so- 

Qj^ called Lifshitz spacetime which has been proposed as the dual to Lifshitz field theories. 

■ For effective mass squared between — (d + z — 1)^/4 and z'^ — {d + z — 1)^/4, we find 

a one-parameter choice of boundary condition type. The bottom end of this range 
corresponds to a Breitenlohner-Freedman bound; below it, the Klein-Gordon operator 
need not be positive, so we cannot make sense of the dynamics. Above the top end 
of the range, only one boundary condition type is available; here we expect compact 
initial data will remain compact in the future. 



1 Introduction 



Recent interest in extending gauge-gravity duality to theories with non-relativistic scahng 
has led to questions about the nature of holography in spacetimes with formally degenerate 
boundaries. Lifshitz- and Schroedinger-dual spacetimes both have this feature, as do the 
warped AdS spacetime solutions in topologically massive gravity. In fact, any spacetime 
which is dual to a non-relativistic field theory must have a degenerate boundary; the non- 
relativistic symmetry means that the time and space components of the metric do not scale 
the same with respect to the radial coordinate. Thus, there will be no single overall conformal 
scaling which we could use to define a nondegenerate boundary metric. 

There are several known ways of dealing with this issue. The most common way is to 
simply define a metric for the field theory at a cutoff surface near the edge of the bulk. This 
method is calculationally practical but does not provide for a theoretical understanding of 
the nature of the degenerate boundary itself. However, as we will show, it is still possible to 
gain further understanding of the available boundary conditions via a cutoff approach. 

In this paper we will concentrate on the Lifshitz spacetime as defined in [1]. We will 
not address the Schroedinger spacetime; for references which do, see [2]-|l]. Hofava-Lifshitz 
gravity (distinct from the Lifshitz spacetime first proposed in [1]) also does not have a 
well-defined conformal boundary; motivated by the preferred time slicing present in Horava- 
Lifshitz gravity, [5] proposed that we define the boundary for all spacetimes by peeling off a 
different factor for each piece of the metric, chosen to ensure that the pieces remain finite as 
the radius approaches infinity. 

In a similar vein, p|l7] define a stress-tensor complex at the boundary of Lifshitz spacetime 
by applying an appropriate conformal factor to each portion of the tensor. Additionally, [7] 
provides a definition of an asymptotically locally Lifshitz spacetime. Unfortunately, in order 
to preserve the asymptotic conditions, it becomes necessary to turn off some possible modes 
when doing holography. |8] sets boundary conditions in a different manner, by first requiring 
that all divergences be cancellable by local counterterms. |9] builds a stress-energy tensor 
at the boundary specifically for spacetimes with z = 2, again imposing a set of boundary 
conditions a priori. [10] studies perturbations in a particular Hamiltonian formulation, again 
imposing boundary conditions and thus limiting the available solutions. 

In this paper, we will explore the possible boundary conditions for a scalar field living 
on a background Lifshitz spacetime, following the approach of [II1[I2]- In [H], the authors 
provide a unique prescription for studying the possible boundary conditions for perturbations 
on a stably causal, static spacetime which are compatible with "good dynamics" . By "good 
dynamics" , they mean it is possible to extended initial data for a wave equation on a static 
time slice to a solution throughout the spacetime, while also maintaining 

• local agreement with wave equation 

• a positive conserved energy 

• time translation and reflection. 
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When they additionally impose a particular convergence condition first defined in [13], [TT] 
shows that their prescription for finding consistent boundary conditions is unique. 

In [12] , the same authors explore this boundary condition prescription for global anti-de 
Sitter space. They find that for a low- mass range just above the Breitenlohner-Freedman 
bound, there is a full one-parameter range of allowed boundary conditions, ranging from 
Dirichlet to Robin to Neumann conditions. These mixed boundary conditions were studied 
in [2] and subsequent papers, and associated with multi-trace perturbations of the dual 



In this paper we will apply the approach of [llj to the Lifshitz spacetime with fiat slicing. 



which is both static and stably causal. We will not propose any specific action for which 
this spacetime is the solution; rather we only consider scalar perturbations living on this 
background. 

We privilege the static time slicing, constructing a spatial boundary for each time slice; 
alternatively, we could imagine the boundary of the spacetime as being at a cutoff surface 
inside spatial infinity. Either approach allows us to consider boundary conditions for a scalar 
field at spatial infinity on each time slice. In the AdS case where z = 1, we reproduce the 
alternate set of boundary conditions as found in [12] . When z > 1, we again find a one- 
parameter family of conditions available in the low-mass range just above a Breitenlohner- 
Freedman-type bound. Above this range, only one type of boundary condition is consistent. 

In Section [21 we review the approach to studying boundary conditions presented in [11] . 
In Section [3l we reproduce the global analysis of the scalar field in AdS done in [12] in the 
Poincare patch. In Section HI we study the boundary conditions for scalar fields in Lifshitz 
spacetimes for rational scaling parameters z > 1. In Section we discuss the comparison of 
our results to previous work, and consider future extensions. 

2 Review of Wald and Ishibashi Procedure 

We will now review the procedure described in [TT]. We are interested in studying the 
behavior of fields solving particular wave equations on a spacetime. If a spacetime is 
globally hyperbolic, then we do not need to impose boundary conditions at the edge of 
spacetime for these fields; just knowing the initial data 00, 0o on one Cauchy surface is 
sufficient to evolve the field throughout the spacetime. 

For spacetimes that are not globally hyperbolic, there is no such Cauchy surface in the 
spacetime. If we want to evolve initial data throughout the spacetime, we must include 
additional information about what can come in through the boundary at spatial infinity or 
at singularities. If the spacetime is stably causal, we can impose boundary conditions that 
produce "good dynamics" . 

We impose stable causality because spacetimes which are not stably causal are arbitrarily 
close to having closed timelike curves. Thus, the perturbations due to backreaction could 
produce causal violations, and prevent us from evolving our initial data in a sensible manner. 



CFT. 



dt'^ dx^ dr"^ 
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Stable causality is equivalent to the existence of a well defined time function which 
increases along all timelike curves in the spacetime. Thus, all stably causal spacetimes have 
well defined time-slicings. If a spacetime is additionally static, then each time-slice will be 
equivalent to any other; we are free to consider the physics on one particular slice S. By 
studying the boundary conditions for fields on a given timeslice S, we thus deduce possible 
boundary data for on the spacetime clS £l whole. 

The procedure developed in [TT] finds the consistent boundary data imposable for a field 
solving a wave equation on a static stably causal spacetime, such that these boundary data 
allow "good dynamics" , as defined in the introduction, for the field in the full spacetime. 
It involves essentially five steps: 

• Pick a static direction t which gives a timelike slicing with slice E. 

• Find the operator A which describes the spatial portion of the wave equation solved 
by (j) and describe its Hilbert space. 

• Delineate the positive self-adjoint extensions Au of A. 

• Use Ajj to define (pt which solves the full wave equation throughout the entire spacetime. 

• Analyze the boundary conditions obeyed by (pt for a given extension Au. 

These steps provide a precise means for extending the operator A to an action Au on the 
spatial boundary of S, which is well-defined since it is a single time-slice. Additionally, they 
allow us to examine the domain of the new Au, and thus the boundary conditions which the 
0f will obey as we propogate the data 0o into future time slices. 

Note that positivity of Au refers to the condition (0|y4[/0) > for all in the Hilbert 
space. Negative norm would indicate an instability in the future evolution of (pt, so we wish 
to avoid it. We now proceed to explore each of these steps in more detail. 

2.1 Finding the operator A 

First we pick out a static direction t. Some spacetimes may allow multiple choices of static 
direction; in the case of AdS, we may choose either global time or Poincare time. We expect 
that boundary condition results should not depend on the choice of slicing. In any case, 
picking a static time direction t privileges a particular time slicing with equivalent slices S. 
We can write the wave equation solved by the scalar field as 



for some operator A that acts only on S. The domain of A is the set of smooth functions of 
compact support which live on E. We treat A as an operator acting on the Hilbert space of 
square integrable functions on E, but with volume measure V~^dT,, where 



= 





(2.2) 



encodes the effect of the time components of the metric. 
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2.2 Finding positive self-adjoint extensions of A 

Given the operator A, we will need to enumerate its positive self-adjoint extensions. We will 
additionally need to insist that A itself to be positive and symmetric; it is this restriction 
which will give us the Breitenlohner-Freedman bound on the mass of a scalar. Here we will 
only summarize the procedure for finding these self-adjoint extensions; in [11], the authors 
explain why such extensions will be related to the available boundary conditions which 
preserve "good dynamics" . 

First, we find the size of the solution spaces /C± whose elements solve 

Aij± = ±^^±, (2.3) 

and are also square integrable under the measure V~^dl^. 

If ]C± are both empty, then there is a unique extension, given by A, the closure of A 
(this is the Friedrichs extension). Thus there is a unique choice of boundary condition in 
this case. 

Conversely, if IC± are the same dimension, then we find the maps U : /C_|_ — t- /C_ such 
that 110+11 = ||f/0+||. Given a choice of U, we then define a particular self adjoint extension 
via 

Au^ = A(l)o + i^+-iU^+. (2.4) 
Note that the domain of Au is larger than that of A: 

Dom{Au) = {00 + 0+ + f/0+|0o e Dom(A), 0+ G /C+} . (2.5) 

Lastly, we check positivity of proposed Au. We have one positive self adjoint extension for 
every U which produces a positive Au. 

2.3 Finding (j)t 

We want to find time-dependent solutions (pt to the equation (12.11) . We can define 

(j)t = cos(A^/'t)0o + A"'/' sin(4/'t)0o, (2.6) 

where 00 and <pQ represent the initial data on our reference time slice So. As shown in [TT], if 
00 and 00 are smooth functions of compact support on Sq, then 0f will be a smooth solution 
to (12.11) throughout the spacetime. 0t matches the initial data on our reference time slice, 
and produces "good dynamics". 

Lastly, we study the boundary conditions satisfied by 0t defined as in (12. 6p . Thus we 
relate the choice of map t/ to a choice of boundary condition on the field 0. 

3 Scalars in AdS 

We now apply the procedure from [11] as reviewed in Section [2] to the case of a scalar field 
in AdS spacetime on the Poincare patch. We will reproduce results in [12], which studied 
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fields on AdS spacetime in global coordinates. [12] did not restrict its results to scalar fields; 
rather, the authors reduced equations for graviton modes and gauge modes to effective scalar 
equations by exploiting the 50(2, d — 2) symmetry. Here, we will concentrate only on scalar 
fields solving the Klein-Gordon equation; we leave analysis of other modes to future work. 

In particular we use the Poincare patch because the Lifshitz spacetime we wish to study 
does not have a good global coordinate system; instead only the analogue of the Poincare 
patch is available. We should also note that AdS global time is not equivalent to AdS 
Poincare time; therefore some of the details of our analysis will differ, but the end result 
describing the boundary conditions at the spatial boundary will be the same. 

Specifically, we work with AdSd+i written as 

dt"^ dx^ dr"^ 



We study the Klein-Gordon equation for a massive scalar in this background, which becomes 
- r^a^V + r'^^^dr [r^~'^dr(t)] + ■ V0 - mg0 = 0. (3.2) 
The Hilbert space consists of functions on spatial slices S which are normalizable under the 



norm 



(02101) = / 0;0iri-^rfrrfx,_i. (3.3) 



s 



Since the wave equation is fully separable, we expand as an integral over plane waves 
labelled by /c; additionally, we will multiply by a power of r designed to remove the single 
r-derivative term in (13.21) . Thus we define ipk{t,r) via 

(j) = r'^ [ dkc^e^'^-'^Mt^r), (3.4) 



where depends only on k and tells us the shape of the wave packet in the x directions. 
Since we are only interested in behavior at the r = boundary, we will only need to consider 
boundary conditions on the ipk- We can now find the operator A by rewriting the Klein- 
Gordon equation (13. 2p as 

dttpk = -A^k, 

- i 

A = -dl + + k^, (3.5) 

where = + {d/2)'^. The Breitenlohner-Freedman bound is saturated for u = 0; A is 
also only positive for masses above this bound. Accordingly we will only study z/ > 0. 
Under the redefinition (13. 4p . the r portion of the norm ( 13. 3 p becomes simply 

{^2\^i)r = j dr^l^tpi^k- (3.6) 

We will now follow the procedure outlined in Section |2l First we will study the eigenvalue 
equation ( 12. 3p with the operator A as in ( 13. 5p under the Hilbert space of functions normal- 
izable under ( 13. 6p . We then use these solutions to construct the self adjoint extensions Au 
of A, and study which boundary conditions these extensions correspond to. 
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3.1 Solving Alp = 

We wish to study solutions of 

All) = AV', (3.7) 
with A as in (13.51) . We can rewrite this equation as 

- r'9,V + - ^) ^ + {''^ - A) = 0- (3-8) 

The solutions to this equation can be written in terms of Hankel functions, which are given 
in terms of the usual Bessel functions as 

Hi^iar) = Ju{ar) + iY^{ar). 

H2,uiar) = J„{ar) - iY„{ar). (3.9) 

Their large r behavior is 



H2,u{ar) = 




(3.10) 



for I arga| < it. The solutions to ( 13. Sp are 

Mr) = Ci^Hi^,{rVX^) + C2v^i/2,.(rv/A^), (3.11) 

where by \/\ — we mean the root with positive imaginary component. This is just a 
parametrization choice and is always available for A = ±2, the case we are interested in. We 
make this choice so Hi will be exponentially damped near r — oo, while H2 will blow up 
there. 

For < 1/ < 1, both y/rYjj and ^JrJy are normalizable near r = under the norm 
(13. 6p . Consequently, as we can see from the large r behavior, the solution Hi will be square 
integrable under the norm (13.61) . The solution 7^2, however, blows up exponentially near 
r — 00 and thus is not square integrable. 

For u > I, neither solution is normalizable, because only y/rJi, is normalizable near r = 0, 
so both solutions Hi and H2 are disallowed. We cannot consider y/rJ^, itself either, because 
it will be a combination of Hi and H2, and will thus blow up near r — )■ 00. 

Thus, for < < 1, we have one (linearly independent) solution to (12. 3p for every value 
of k and each of A = ±i, and thus the dimension of both ]C± is 1. For u > 1, IC± are both 
empty. These results match with those of [T2] . 

3.2 Finding the extensions Ajj 

We now use the results of the previous section to build the self-adjoint extensions Au. The 
arguments follow similarly to those in [T2] . 
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First, for the case > 1, the only available extension of A is just the closure of A, A. 
That is, for u > 1, we simply extend the action of A on each slice S to the spatial boundary 
of S at r = 0. 

For perturbations whose effective u satisfies < < 1, both and /C_ are one dimen- 
sional. Specifically they are spanned by 'ip±, defined from (13.111) as 

V'i = i^x=±i = C^V^H^,, (rV±i - , (3.12) 

where — W refers to the square root with positive imaginary component, and the C± = 
C are positive real constants chosen such that \\ip±\\ = 1 under the norm (13.61) . It can be 
shown that C is the same for both ip^ and ip^. 

As the sets /C± are one dimensional, we can write all maps which take /C+ /C_ in the 
form 

[7V'+ = e^>_, (3.13) 

for —IT < a < n. In other words, we can only introduce a phase into the map, but the 
phase can take any value. Note that the transformation U takes the element ip^ of /C+ and 
associates it with a particular element, e*"'0_, in the space /C_. Our only definition of the 
action of U is as given in (I3.13p . 

For a particular choice a, the map corresponds to the extension of A whose operation on 
ip is given by 

A^ijj = Aipo + i^p^-ie^'^'ilj.. (3.14) 



3.3 Interpreting the boundary conditions 

Again following [12], we will need to understand the asymptotic behavior of 



(3.15) 



near the boundary r = 0, for < < 1. Expanding ipa near r = for u not a half integer, 
we find 



_|_ J, 



1/2+u 



1 + e*" 
1 



V-i -ky 

i cos(7rz/)r(— z/) 

TT 



(3.16) 



+ 0(r^/2-.)_ 



We see from this form that for l/2<z/<l, we allow for solutions which blow up as r — t- 0. 
They still obey prescribed boundary conditions, however; specifically, solutions for which 
the r^/'^~'' coefficient is nonzero correspond to generalized Neumann or mixed boundary 



7 



conditions, as described in [12]. In order to compare to their analysis, we now rewrite the 
coefficients in square brackets in (13.161) as 



Recalhng that the chosen square roots are those with positive imaginary components, we 
can reduce these coefficients to 

au = I + exp — iu arctan — ^ , (3^ = 1 + exp (^a + iu arctan — ^ . (3.18) 

Both ay and [3^ range from to 2. Consequently their ratio, Pu/a^, ranges from to oo. 
We can thus safely label our choice of self-adjoint extension by the this ratio, rather than 
the value of a. One reason for this choice of labelhng here is that that the particular value 
of a which produces a given ratio f^vjoiy depends on k] this is an artifact of the Poincare 
coordinates we have chosen to work in. 

More importantly, the ratio fiyjoiy corresponds most directly to the boundary conditions 
implied by the choice of extension. As in [12], further analysis is simplest for the case 
z/ = 1/2. We must use a different expansion for the Hankel functions valid for z/ = 1/2, but 
the analysis follows exactly as in [T^ . 

In this case we find that Pv/ay controls the ratio drip/ip at the boundary r = 0. Thus 
ttu = 0, or f3u/aiy = oo, corresponds to regular Dirichlet conditions, P^/ay = corresponds 
to Neumann conditions, and other values correspond to mixed Robin conditions. 

For more general values of < z^ < 1, [12] call the cases Pu/a^ = oo, generalized Dirich- 
let, Neumann conditions, and all other choices are termed generalized Robin conditions. 

The domain of a particular extension, as defined in (12. 5p . does determine the boundary 
conditions for the time-propagated solution (pf That is, even though we begin with compact 
data ipo, ipo, the definition of 0t in (12. 6 p will ensure in general that future slices have instead 
the boundary conditions of the full domain (12. 5p . 

Recalling the definition of ip in (13. 4p . we can recover the boundary conditions for (pf We 
find the behavior 

where again Pu/a^ controls the boundary conditions, and # represent numbers not dependent 
on the choice of a. We note that not all boundary conditions may be possible for a given 
k, u] we have not tested for positivity of the extension. In fact, [12] found that only a partial 
range of choices for Py/a^ produce a positive extension. We leave out these details here; the 
interested reader can see [12] for the precise ranges involved. 



4 Scalars in Lifshitz 

We now consider the Lifshitz spacetime in c? -|- 1 dimensions, with metric given as: 

dt'^ dx^ dP 
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where x runs over d — 1 dimensions. This spacetime is static and also stably causal. The 
boundary in these coordinates is at r = 0. We will not discuss the tidal singularity present 
at the center of the spacetime where r = oo. Indeed, we will see later in equation fl4.8p that 
the natural norm on scalar fields will require that these fields be strongly suppressed near 
the tidal singularity. 

For z > 1, multiplying by a conformal factor r^^ produces a degenerate metric at r = 0. 
As discussed in the introduction, this behavior is by design, since this spacetime is designed 
to be dual to a field theory with different scalings for time and space. For our purposes, we 
will privilege the time slices S defined by constant t, which are d dimensional and have a 
well defined non-degenerate d — 1 dimensional boundary at r = 0. 

Before we analyze the operator A as defined in (12. ip for the Lifshitz spacetime, we wish 
to highlight the issue of the Hilbert space in question. As mentioned above, the operator A 
acts on the Hilbert space of functions which live on E and are square integrable with respect 
to the measure V~^dT,. Specifically for the Lifshitz metric (14. ip . we find V = r~^, and thus 
the appropriate measure is 

y-^dS = r^dS = r^~^ dr dxd-i (4.2) 

where dxd~i is the volume form for the slice in Lifshitz geometry. The inner product 
is given by 

(02|0i)= / (t>*2(i>ir'-''drdxd-i. (4.3) 

We can see that the set of functions included in the Hilbert space does depend on z. 
The Klein-Gordon equation V^V^0 — mQ</) = for the metric (14. ip is 

- r'^'d1(j) + r'^+'dr [r^-^-^a,^] + r^V • V0 - m^^ = 0, (4.4) 

where V is the gradient on the R"'^-'^ slices. We expand cf) as 

(j) = rP J dkc-^e^^-^'Mt^r) (4.5) 
d + z-2 

where only depends on the value of k. The Vfc only depend on the magnitude of the 
momentum k, as we can infer from the form of (14.41) . With this expansion we can now 
rewrite (14. 4p in the form prescribed in (12. ip : 

d^Mt,r) = -AMt,r), (4.6) 

where the operator A is given by 

A = -r2(i-^) 



(4.7) 
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We now would like to write the norm (14. 3p in terms of the ipk- The e^^'^ are orthogonal but 
not normalizable. Since we are interested in the behavior at the boundary r = 0, we will 
only study the behavior of the i/jk{t,r) pieces, assuming that we can arrange wave packets 
which are normalizable in the x directions. In terms of tpk{t,r), the relevant portion of the 
norm (14. 3 p becomes 

(02|0i). = J drr'^-'r2,kit,r)^,,,{t,r). (4.8) 

We now have set up our problem: we want to study boundary conditions at r = for 
ijjk{t,r) which solve (14. 6 P with A as in (14. 7p and are finite under the norm (14. 8p . In the 
remainder of this section, we will study the behavior of solutions to the complex eigenvalue 
equation (12. Sp . then define the self adjoint extensions Ajj of this A, and finally study which 
of these extensions are finite and what boundary conditions they correspond to. 



4.1 Finding the solution spaces K.± 

We now return to the equation (12. 3 p with operator A as in (14. 7p in order to find the solution 
spaces /C±. Since we only want to understand which solutions to (12. 3 p are normalizable 
under (14. 8p . we actually do not need to solve for the full eigenfunctions ip±] instead we will 
only need to understand their behavior near the boundary. For ease of notation, we will 
actually analyze the behavior of the eigenfunctions for which Aip = Xijj for generic rational 
values of z, and general complex eigenvalues A. We also drop the subscript k. We rewrite 
this eigenequation as 

r^d^ip - (''^-\)'^- ^^^V + Ar^'V = 0, (4.9) 



u^^ml+i^±4-^] . (4.10) 



where u is given by 



Note that 1/^ = corresponds to the equivalent of the Breitenlohner-Freedman bound in this 
space, mQ = —{{d + z — l)/2)^. Considering only z/^ > ensures the operator A is positive 
on our Hilbert spacej^ Consequently, we will only consider u > in our subsequent analysis. 

Equation (14. 9 p is a second order linear differential equation with two singular points, 
r = and r = oo. As such, there are two linearly independent solutions for any given A. We 
are only interested in solutions which are square integrable under the measure (14. 8p : thus 
we will examine the behavior of solutions near each of the singular points for such square 
integrability. 



^Usually fields below the Breitenlohner-Freedman bound arc associated with instabilities; |T5] indicates 
that the situation is different for Schroedinger spacetimes. Lifshitz apparently is more similar to AdS; we 
think the positivity bound on A for here indicates that massive scalars below the BF bound in Lifshitz 
will suffer from an instability. 
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The singular point at r = is regular as long as 2 > 0, which includes the range we 
study: z > 1. Considering a power series solution ip = '^^anr'''^"' to (14.91) . the equation 
becomes 



=ao 



b{b 



1 



r + ai 



+E 

n=2 

00 

n=2z 



an{b + n){b + n 



{b+l)b-u- 



1, 



1 



3+1 



1) - (z/^ - -)a„ - k'^an-2 



(4.11) 



a'n-2zr 



b+n'] 



for half- integer z. The coefficient for each power of r must vanish separately, so the indicial 
equation becomes 

b{b-l) - v"^ + - =0. (4.12) 



For rational z which are not half integers, we may derive the same effective indicial equation 
by changing variables to p = r^/"^ for m such that mz is an integer. This equation gives the 
leading behaviors of the two solutions for %p near r = as 



(4.13) 



assuming u is not itself zero or a half integer, u must be nonnegative in order to assure 
positivity of the operator A on the desired Hilbert space. 

When u is zero, the two roots of the indicial equation (I4.12p coincide, and thus the 
leading behaviors become r^/^ and r^/^ logr. For u some nonzero half integer, the two roots 
will differ by the integer 2z^. Here, the leading behavior as r approaches zero will be r^/'^-'^^ 
but there may be logarithmic terms starting with r^/'^+'^ logr in the full expansion. 

A solution ilj{r) will be square integrable near zero with measure as in (14.81) when the 



leading behavior of ■^^2^,22-2 jg a > —1. The solution ip^ has behavior near r = such 

that 

= ^2z-2+2i/+i. square integrable if > -z. (4.14) 

As we are working only with u > 0, this solution will always be square integrable near zero. 
For the solution when u is strictly greater than zero, we find 



„22-2-2i/+l. 



square integrable if 



u > —z. 



(4.15) 



as the possible logarithmic behavior is subleading. For the case z/ = 0, we wish to study the 
leading logarithmic behavior 

^2.-2^2(1/2) ^2^ (4.16) 



which is integrable near zero for any z > 1. 
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These results indicate that for < u < z, both solutions i/j^^^ are square integrable 
near zero. For u > z, only the solution remains in the Hilbert space, regardless of the 
eigenvalue A under consideration. 

In order to find the possible square integrable solutions to the eigenvalue equation (14. 9p . 
we must also consider the behavior of these solutions near the irregular singular point at 
infinity. Again, as near r = 0, we expect two solutions each with different behavior. 

The point at r = oo is an irregular singular point with rank z. Thus there must exist 
formal solutions, or asymptotic series, of the form 



ip = exp 



n=l 



Q{r), Q{r) = 



air 



b-l 



(4.17) 



1=0 



For more details on such solutions, see e.g. [16]. Plugging this form into (14.91) and dividing 

by exp [Y,n=i '^nr'^] gives 



=dlQ + 2^nCnr''-^drQ 



(4.18) 



n=l 



+ 



^^n(n — l)c„r'^ ^ + j ^^nc„r 

n=l \n=l 



Q. 



We will start by setting to zero the coefficient of the highest power of r, as these are the most 
dominant terms near r = oo. The largest power in Q is r*; for drQ the largest power is r^~^, 
and for d^Q it is r''"^. Thus the largest overall power present in (I4.18P is r^^"^"*"*, arising 
from the A term and the highest term in the squared sum in the second line. Specifically 
this contribution is 

r2--2+^ (^2^2 ^ A) = 0, (4.19) 

and it must vanish independently. Note that if 2; = 1, the —k"^ term is of the same order 
and would contribute; this is precisely the AdS case studied above. Here we will restrict our 
attention to ^ > 1. As for the behavior near r = for rational z, we can find the behavior 
for rational, noninteger z here by changing variables to p = r^/*" for m such that mz is 
an integer. We again recover equation (14.191) for the coefficient of the leading power in the 
exponent of ip. 

The behavior of ip near infinity is controlled by Cz, which must be a root of equation 
(I4.19p . Labelling the two roots of this equation c^^-t, we can write the two allowed behaviors 
of ijj near infinity as 



i/ji oc exp [cz,+r' + 0{r'-^)] (/ + 0{r^~^)) 
1P2 oc exp [cz,-r' + 0{r'^^)] (r^ + C(r^^i)) . 



(4.20) 



If the real part of is positive, the solution blows up exponentially and will not be square 
integrable under the measure (14.21) . Conversely, if the real part of Cz^± is negative, then the 
solution is exponentially damped and will be square integrable. 
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Since we are interested in the solution space for (12. 3p . we now examine the roots of fl4.19p 
for X = ±i. We find 



In either case, the solution 1/J2 built on Cz,- is square integrable, and ipi built on Cz,+ is not, 
for every rational z > 1. 

We are searching for solutions to fl2.3p which are square integrable under measure (14.21) 
over their full range. The only candidate solution near r = 00 is ?/'2- This solution will be a 
linear combination of those near r = 0, denoted 'i/'t,i in (14.131) : 



for constants Cf, C|. For < u < z, both tp-^ and ipi are appropriately square integrable 
near r = 0, so in this range the solution 1^2 itself will be square integrable with measure (14. 2 p 
from < r < 00. 

Conversely, for u > z, the solution ip^ grows too fast near r = 0. Thus ip2 will not be a 
valid square integrable solution to (12. 3p unless Ci = 0, which will not occur generically. 

In summary, for < u < z, one generating solution to (12. 3p for each of ±i is square 
integrable. Since we can multiply this solution by any phase, we have a one- dimensional 
set of solutions for each of ±i, and thus we will find a one-dimensional alternative set of 
boundary conditions available for modes with this range of u. For u > z, no such solution is 
available, and only one type of boundary condition is sensible. 

4.2 Boundary conditions for i/j when < u < z 

We have shown that there is a one-dimensional space of possible boundary conditions for 
scalar operators whose mass satisfies u < z. We now wish to explore what boundary con- 
ditions these are. The added complication of the factor z in the equation (14. 9 p will prevent 
us from finding as explicit a form for the boundary conditions as we could in the AdS case. 
In particular, we will not be able to find the coefficients Cf, C| which determine the precise 
behavior of the function ip2 near the boundary r = 0. However, although we will not be able 
to find these values, we can still discuss the available behaviors near the boundary. 

As in the AdS case, the maps U from the spaces /C± for square integrable solutions 
X = ±i are of the form (I3.13P and thus just introduce a phase. Consequently we again define 

V^Q = i^2,\=i + e"^tp2,\=-i- 

In (I4.13p . we find the behaviors of ip^^^ near r = to be 1/2 ± i^. For u not a half integer, 
we can use these behaviors to write the behavior of ipa as 




(4.21) 



^2 = CtV't + C'iipi, 



(4.22) 



00 




(4.23) 



n=0 
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where the coefficients c„^q, bn,a depend on the choice of map, and control the choice of bound- 
ary conditions. Unfortunately our inability to match the solutions ip^^i near the boundary 
r = to the solutions ipi^2 near the tidal singularity r = oo means we cannot find the 
ranges of the coefficients Cn,o, bn,a- Similarly we cannot test within this range to find which 
boundary conditions correspond to positive choices of extension. 

However, we still expect that the choice of extension, or of a, corresponds to a choice for 
the ratio CQ a/bo,a- We still define (pt via (12. 6 p and expect its boundary conditions to explore 
the full domain of the choice of extension determined by Co,o/feo,o- 

Using the relationship between and ip as in (14. 5 p , we find the behavior of (pt near r = 
to be 



The leading term here goes like 7--'^+("'+2-i)/2 As choices of boundary conditions are available 
for scalar fields whose u satisfies < u < z, we find that the exponent of this term will be 
bounded below by [d — z — l)/2. Thus, for Lifshitz spacetimes in ci + 1 dimensions where 
z + 2 > d + 1, one can choose that some scalar field will have a blow-up near the boundary 
r = 0. 

This behavior could occur when the mass is as far above the Breitenlohner-Freedman 
bound as we can achieve and still have a choice of boundary condition behavior. We found 
that a choice of boundary conditions is possible whenever < u < z, oi when ttlq satisfies 



We expect that scalar fields with values of near the top end of this range will have a 
boundary condition choice resulting in growth near r = 0, as long as the Lifshitz spacetime 
under consideration has z + 2 > d + 1. In fact, spacetimes above this limit have been 
considered previously but been found to have other sicknesses; see [7j and references for 
details. 

The second series in (pt controlled by co.q, begins down by a power r^". Again, as we can 
pick alternate boundary conditions for values of u up to z, we can see that these conditions 
correspond to a further generalization of the Neumann, Dirichlet and Robin conditions. In 
particular in the case u = n + 1/2, we might guess that the alternate boundary conditions 
here relate to the derivative dl!"ip. 

5 Summary and Future Directions 

We have studied the boundary conditions available for Lifshitz spacetimes which are solutions 
to Einstein gravity plus unspecified matter. We find a one parameter family of boundary 
condition types for scalar fields with effective masses satisfying (14.251) . In the AdS case, we 
label these conditions by the ratio a^/ Pv, where multiplies the leading behavior, and pi^, 
the subleading behavior. For the Lifshitz case, we define similarly 6o,a; co,a- A choice of ratio 



oo 




(4.24) 




(4.25) 
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du/ f^u, or &o,q/co,q, corresponds to a choice of self-adjoint extension or boundary condition 
type. These types range from a generahzed Dirichlet condition, through a family of mixed 
conditions, to a generalized Neumann condition. We summarize our results in Table [TJ 





1^ = 


Q < U < Z 


U > Z 


mass 








(pt leading order 


iy,{d-\-z l)/2 I^Qg ^ 


^-v+(d+z-l)/2 




(f)t subleading 


^{d+z~l)/2 


^u+{d+z-l}/2 


r^V+[d+Z~l)/2 



Table 1: Summary of results: the leading order behavior of (pt is controlled by Uy for AdS, 
or 6o,Q for the Lifshitz case. The subleading is controlled by [iy or cq^q. Thus the boundary 
conditions are set by the ratio a^jfiy^ or ho^a/co^a- For z/ > 2;, no choice is available. 



Additionally, following [12], we can define a conserved, positive energy for each choice of 
boundary condition, parameterized here by — tt < a < tt: 

E{ij) = J drr^'-^*ij + J drr^'-^r (^^ + #2,A=i - ^e'>2,A=-i) , (5.1) 

where A is given in fl4.7p . and we have used the definitions in 03.141] with ip2,x as in fl4.2UI) . 

Below the mass range fl4.25p the Laplacian operator is not positive; the lower end corre- 
sponds to a Breitenlohner-Freedman bound. Above this mass range, we cannot choose the 
type of boundary condition; only a Dirichlet-type condition is available. Importantly, we are 
able to reproduce the results for scalar fields found in [12] for the z = 1, that is the AdS, 
case, albeit in Poincare coordinates. 

The authors of [12] also found alternate boundary condition options for gravity and 
gauge fields. The authors accomplished this by expanding these fields in modes in global 
coordinates, resulting in a set of effective scalar equations with an extra parameter a. As 
Lifshitz spacetime does not have a good set of global coordinates, we have not done a similar 
expansion. 

We have also not investigated the effect of the tidal singularity (at r = 00 in our coor- 
dinates). The method of exploring boundary conditions proposed in [TT] should work for 
curvature singularities as well as for boundaries at infinity; it would be interesting to study 
the method for tidal singularities as well. Although we have not done this analysis, we 
already see that the solutions to (12. 3p which are square integrable are in fact exponentially 
damped as exp(— r^) in the region r — )■ 00. 

It would also be interesting to compare the alternate boundary conditions we have found 
here to previous work on holography in Lifshitz spacetimes. Particularly, we wonder about 
the relation to a certain mixed boundary condition required in [7]. Additionally we have only 
considered Lifshitz spacetimes as solutions of Einstein gravity; interesting recent work [17] 
suggests that Hofava-Lifshitz gravity could be a natural ground for holographic duals to 
Lifshitz field theories. 

One can also ask what these alternate boundary conditions correspond to in the putative 
dual Lifshitz field theory. In the AdS case, a choice of mixed boundary condition corresponds 
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to adding a multitrace operator in the field theory, as studied in "designer gravity" as in [H] . 
It would be interesting to develop a similar "designer Lifshitz gravity" story. 
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